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TECHNICAL MEMORANDUM No. 602 

WRINKLING OF REINFORCED PLATES SUBJECTED TO SHEAR STRESSES* 

By Edgar S§ydel 

Resume 

It is common practice to stiffen flat plates toy ribs or 
special shaping (corrugations). Now in many load cases the 
principal stress in a plate is the shear which, upon reaching 
a certain limit, the so-called critical load, forces the plate 
to wrinkle. 

When these reinforcements are appropriately applied the 
stiffened plate can be treated as a homogeneous plate (orthogo- 
nal-ani so tropic) with uniformly and closely spaced reinforce- 
ments (different in both directions), and the formula for buck- 
ling in shear is deduced from a differential equation of the 
fourth order. By introducing the characteristic value of the 
orthogonal-anisotropic plate and applying further appropriate 
parameters the solution can be carried out for a plate strip of 
great length by arbitrary plate stiffeners. The edges may be 
freely supported or rigidly constrained. 

*"Beitrag zur Frage des Ausbeulens von versteiften Flatten bei 
Schubbeanspruchung. " From Jahrbuch 1930 der Deutschen Versuchs- 
anstalt fur Luftfahrt e.V., Berlin-Adlershof - Verlag von R. 
Oldenbourg, Munch en-Berlin 1930, pp. 235-254. 
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, I; n t ; ;r :; p-;d, u -c ."ir-i o n 

•. The buckling strength in shear i and the vibration) of long, 
thin, isotropic plates supported; at the long edges has been 
treated quite exhaustively by R. V. Southwell and S. W. Skan 
(Reference l). St. Bergmann and H. Reissner (Reference 2) ex- 
tended Southwell's theory to include the case of anisotropic 
( or thogonal-ani so tropic) plates with vanishing and comparably 
low bending stiffness longitudinally. Their data are included 
in this report. 0. Schmieden (Reference 3) treats the case of 
equal longitudinal stiffness in bending and distortion, but his 
methods were incomprehensible to me and his data are confirmed 
only for the case of infinitely low distortion stiffness. 

M . T. Huber (Reference 4) established the theory and the 
general differential equation" of bending in or thogonal-ani so- 
^tropic plates. 

In the present report the problem of buckling strength and 
lobed form of failure of anisotropic plate strips is carried to 
a certain conclusion and the data for arbitrary stiffness in 
bending in two directions at right angles to each other, and for 
arbitrary distortion stiffness is exhibited diagrammatic ally. 
For the suggestions on this report and helpful advice I am in- 
debted to Professor H.. Reissner. 
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problem 

1. The system (the reinforced plate) .- The thin plate (met- 
al or plywood) is a structural element which finds extensive use 
in a ir plane construction. Among others it is intended to form 
a substitute for "bracing and for that reason it is in many load 
cases chiefly subjected to shear. The low stiffness of the thin 
plate which induces bulging (lobing, wrinkling) even under rela- 
tively low shear stresses, and the inability to absorb compres- 
sion and bending stresses, necessitates reinforcements. This 
stiffening can be accomplished with riveted. or .glued strips 
(ribs), or by special shaping (corrugated plate), usually in 
one or two. directions. 

Wrinkling is, an 'elastic form change in which deflections 
and distortions occur in the originally -flat, plate center, and 
depends therefore on the stiffness in bending and distortion of 
the plate. . 

This stiffness (with respect, to unit length^) of the homoge- 
neous isotropic plate, • •being the same at every point and in ev- 
ery direction, is higher at the reinforcements than at the non- 
reinforoed points. In addition, the ■ bending .stiffness (depend-, 
ing on shape of. reinforcement ) is usually different in differ- 
ent directions of the stiffened plate. 

If the reinforcements are equally spaced and under certain 
circumstances of 'equal thickness, the stiffness (with respect 
to unit length) is periodically variable along the two axes. 
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A flat plate may "be visualised as taking the place of a corru- 
gated sheet (located in the ideal median plane of the corru- 
gated sheet), whose stiffness corresponds at every point with 
that of the corrugated sheet. In this case the stiffness again 
is subject to periodic changes - here, continuous - while in 
the plate with reinforcing strips the cross sectional inertia 
moments and through them, the stiffness changes irregularly 
(at least, when making the usual simplifying assumptions regard- 
ing the stress of a beam in bending). Thus, the stiffened plate 
represents, even if it consists of isotropic material, a non- 
homogeneous system, "because the elastic properties of the plate 
are not the same at every point. But we speak of a mean bend- 
ing stiffness in one and in the other direction (parallel and 
at right angles to the direction of the reinforcing strip). 
Hence, we must distinguish between the mean bending stiffness 
B x , and the mean bending stiffness D 2 , " with respect to two 
axes given by the direction of the reinforcements, and which 
may be assumed as the x and y axes of a system placed in the 
(ideal) median plane of the reinforced plate. Hereby D x is 
the bending stiffness corresponding to the stress in bending 
o* x in the x-direction, and D2 conformal to the normal ten- 
sion c In this manner we arrive at an idealized homogeneous 
plate system whose elastic' properties are different in two di- 
rections, which, in fact, is conformal to a flat nonreinf orced 
plate of orthogbnal-ani'sotropic 5 material (originally, isotropic - 
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or, better, quasi-i so tropic material when rolled, assumes suoh 
qualities to a slight extent). 

• Now, when the periods of the stiffness-change of an actu- •• 
ally reinforced .plate are sufficiently small, the elastic sur- 
face of this nonhomogeneous plate in bending (and distortion) 
deflection differs only slightly from the elastic surface of the 
corresponding homogeneous orthogonal-anisotropic plate, and we 
may speak of a quasi -homogeneous system. Hereby it should be 
borne in mind that the periods of stiffness-change must not only 
be small with respect to the plate length, but also be so small 
that- the elastic, surface, of the ideal system conforms -'to- the' , 
actually existing . system. " ■■■ 

Disregarding the transverse elongation with its trifling 
and frequently negligible effect, the ! bending stiffness is the 
product of the Young's modulus of the material and the (mean) 
inertai moment of the reinforced plate section per unit length. 
The determination of the distortion stiffness requires an exam- 
ination of the elastic form changes of the reinforced plate, 
where it is assumed that only shear moments (no bending moments) 
act in every intersection parallel to the x-axis and parallel 
to y-axis, due to the load being applied to the edges. As a 
rule, the stiff eners do not materially raise the distortion 
stiffness of a flat plate when the reinforcing strips are small 
in comparison to their spacing,- although the edge attachment 
of a corrugated plate seems to affect its distortion stiffness 
considerably. 
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The present report is confined to the rectangular plate of 
(infinitely) great length (in the x-direction) and of 2a width 
(in the y-direction) . The stiffeners are parallel to the edges, 
and the plate thickness 6 is small with respect to the width. 
At the edges y = ± a the plate is freely supported or rigidly 
constrained. 

2. .Loading .- Such a strip is subjected to an evenly distrib- 
uted constant shear t , per unit length, which, acts along the 
four edges in the median plane of the plate and is in outside 
equilibrium. The result is a constant shear stress in the whole 
quasi-homogeneous plate. 

While t is low the plate remains flat, but as the shear 
increases and finally becomes maximum, the plate reaches its so- 
called shear strength t^ r , and upon exceeding this critical 
limit, the plate wrinkles. Figure 1 shows the basic form (con- 
tours and cross sections) which the elastic surface assumes by 
incipient wrinkling. (This elastic surface is discussed in a.' 
subsequent section (III) . It depends on the ratio of bending - 
and distortion stiffness and the boundary conditions.) The 
stress attitude of the plate changes, bending and normal stresses 
occur, the plate attempts to assume a new attitude of equilib- 
rium, but fails if the ultimate strength is exceeded or local 
buckling sets in. We begin with the determination of the shear 
strength tjj. r . 
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II. T -e s t Data 

5. The shear strength. -.- Before going into details, we wish 
to state the result. 

The sheaf strength t kr is computed from one of two formu- 
las': 



W ~ c a a3 



^• Pl 8 , q>3 " (la) 



t kr = c b — . (lb) 



a" 



or, according to Bergmann and Reissner 



a 2 



ot/V "" '• (l»a) 



with . D x 

- P 



> (2) 



Here a = half width of plate strip, 

D 1 = mean stiffness in bending in direction of 
x (length), 

D 3 = mean stiffness in bending in direction of 
y (width) , 

D 3 = mean stiffness in distortion (when disregarding 
the transverse .elongation) f 

of the reinforced plate, the stiffness is expressed in units of 

cross s e ction al wi dth,, that, is, in kg cm, while c a and c^ are 
*See footnote, next page. 
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absolute coefficients contingent upon 

A _ yp, p s ,,yr ( ) 

D 3 ' M- 

and which become apparent from Figures 2a and 2b. . If £ > 1, 
these figures show that it is preferable to use equation (la) 
or (l ! a), but if £ < 1, that formula (lb) or (l'b) would be 
more suitable. When the bending stiffness vanishes in the 
x-direction, that is, if T) 1 = 0, then S- = 0, and formula 
(la) yields t^ r = » x 0, or an indeterminate value, but (lb) 
and (l'b) give in this case 



and 



D 2 D 3 D a y M- 

t krf = 11.71 ■ ■ ^ = 11.71 a 3 ■ by free support 



tj cre = 18.59 fa" ~ 18 - 59 — b y restraint 



on the edges y = ± a. When the bending stiffness disappears 



*( Footnote from page 7) 

Considering the idealized Poisson's ratio, u x and u y in the 
x~ and y-direction (compare section VI, formulas (6) and (7) ) 
would yield the more accurate values for the bending stiffness 
in unit length (EJ) T 

D, = in x-direction, 

D = C s J)y . in y-direction, 
2 l- u x u y 

while D 3 is given by 

? n - v (EJ) y + u (EJ) * + 4 f PrT ) 

where 2 (GJ) X y is the distortion stiffness of the orthogonal- 
anisotropic plate with respect to the cross-sectional width 
"one." The definition of ^ x and is given in section VI. 
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in the x-direction the shear strength t^* is therefore not 
zero (but D 2 = 0 yields t^ = 0). On the other hand, if the 
distortion stiffness is D 3 -= 0, then S- ■= oo an d, according 
to (la) and (l'a) 

t krf = 8.125 djh^k- = 8.135 ~~^> by f ree support, 

and 4 . ^ 

.t kre = 15.065 - — 1 2 3 = 15.065 -2-g — by rigid restraint, 

while (lb) for the time being yields an indeterminate value 
c° x 0- 

If value $ is very large or very small, or, at any rate, 
much higher or lower than "one," one may figure with the c a 
value for S- - » an d c b for # = 0, where the t^ figures 
will be somewhat too low, the error committed here can be seen 
on Figure 3a where the quotient Ca/cg^^o.) is plotted against 
The shear strength . (for ■& > 1) is computed as 



t kr = -~- T 8.125 . ■ ■ _ . by free support 

°a& = ») a 2 

and 4/T "7" 

Co 7 A D 2 

= — 15.065 ■ . 1 - by rigid restraint 

c a(£ = c*) a 2 

where c a /c a (# =O0 ) slightly exceeds "one," and on Figure 3b with 

•Bergmarm and Reiesner first resolved the problem by assuming 
D x =0. Their value — . — 

(JJ>, - 8.285 V 2 X_yQj, 

■ £> 2 /'min ' a2 

is in exact agreement with the above . formula. In their second 
report they showed the initial tangent of the function 

Ob = °b ( *" 1 - 
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°b/ c b(£=0) plotted against l/& . Thus the shear strength for 
£ < 1 becomes: 



t k = — -5h , 11.71 . V 2 ^ (by free support) 



and = .._f.lL.,_ 18.59 ^JkJk (by rigid restraint) 

c b(£=0) . a 

quotient c- D /c-k ($ = o ) being slightljr higher than "one," when ti- 
ls materially lower than "one." When & = 1, so that p = u 3 * 

*The case = 1 was examined by Southwell, and our data, al- 
though established for the isotropic plate, apply for certain 
oases of- the anisotropic plate as well. Schmieden first treats 
the case p = m(=A) and when resolving the problem, confines .. 
himself to very small P (=A) values. He assumes the k value, 
which he denotes by X, and which is inversely proportional to 
the wave length of the elastic surfs.ce, to be inversely propor- 
tional to t/P (= t/A~] in the case of critical shear load, so 

that Kfcr yp~ (= ^vS.) is constant, (i.e., unaffected by P) . From 
the specie! case p = \x < 1 he proceeds to the more general case 
|jl = 0 , p ± 0 and, assuming vanishing plate thickness, publishes 
the result without giving a correct derivation. When P = M- = A, 

then i9- = 1/ J~k\ when A 'moves toward zero , £ goes toward in- 
finity. The term * equals a J £ a = "/O'a* and - e Q. uals 
TT/c a ' in the shear strength case; here c a ' .depends on & only. 

Thus it may be proved that K jCr a/p~ is actually constant for a 
certain £ value (say, for ■&■ =»). With Schmieden's assumption 
p = u, = oo 5 of course, is synonymous with p = |i = 0; that is, 
with the supposition that bending stiffness D a and distortion 
stiffness ' D 3 disappear, in this case the shear stress is alto- 
gether unstable. To confine the computation of c a in the shear 
strength formula to this case alone would be superfluous, be- 
cause the formula yields t^ = 0 anyway. The proof that fac- 
tor c a (?/hich is affected by arrived at for' $ = <», can be 
introduced for high, but finite & values, becomes apparent in 
this report. According to Figure 3a, the error committed (say, 
for & = 20), do'e.s not exceed 5. 6$ when substituting factor o a 
for & = oo (instead of & - 20), while with & = 1, the discrep- 
ancy is nearly 50$. When the reinforcements are arranged in 
length and width as Schmieden assumes (u = 0) , the supposition 
of vanishing plate thickness equivalent to that of vanishing dis- 
tort:, on stiffness, corresponds to the case S- - co, which (with 
u = 0) yields a shear strength different from zero. Schmieden's 
computed numerical values agree with our data for & = °°. 
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then 

°a = c b- 

To make a stable shear stress at all possible requires a 
certain bending stiffness D 2 (>0) although the shear stress by- 
vanishing Dj_ or Eb is stable at the beginning (as long as 
"k < ^kr ) • 

To enable us to make comparison with Roissncr*s and Berg- 
mann's latest experiments in which they prove that c^, when 
approaching d = 0, is a regular function of £ 2 = p/m 2 and 
c a when near l/& = 0 a regular function of l/# = \^/Jp7. 
vtq plotted in Figures 4b and 5b," the c a factors against # 2 . 
(The c-k factor is the same as in Figure 2b, the only difference 
being the abscissa. )' The dotted 'secants through the points with 
abscissas # 3 = 0 and # 3 = 0.04, in Figure 5b, denote the di- 
rections of the tangents in the points with the abscissa $ z = 0. 
Figure 4b and 5b are alike but for the larger scale of # 2 and 
its range confined to 0 = & 3 = 1. 

In similar manner Figure 4a shows the ■ c factor plotted 
against l/# (like in Figure 2a but for the abscissa), and 
Figure 5a, c a 'plotted against l/£ for 0 = l/# =1. To ap- 
proximate the tangent direction in the points with abscissa 
l/d- = 0, the dotted secants are defined by the points with ab- 
scissa 1/^ = 0 and l/» = 0.025 (* = » and & =40). (The 
abscissa scale is different from Figure 4a. ) 

On comparison, we now find that curve C'b^ 2 ''' in fa- 0 * 
from P = 0 (Dj. = 0) to about £ 2 = p /V s = 0.25, that is, 
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& =0.5, can "be substituted by its initial tangent, and that 
curve c a (l/#) from u = 0 (D 3 = 0) to approximately l/S- = 
M7*/p = 1 can be substituted by' its initial tangent. As a re- 
sult eft, in the neighborhood of P = 0, can be developed in 
power series according to & 2 , and c a , on approaching u = 0 
in a power series according to l/£, as proved analytically by 
Bergmann and Reissner. 

The diagrams can be summarised as follows: 

Introducing the constants 



°ao 


8.125 


and 


15.065 


°ai - 


5.64 


it 


7.685 


c a2 = 


0.6 


1! 


0.6 


°bo = 


11.71 


II 


18.59 


°bi = 


2.155 


II 


5.86 


C b2 = 


0.7 


II 


2.3 



we have for pin- jointed and constrained support 
t VT , a 2 



lei /-» i _ 1 . _ 1 

i7= - G a _ G ao " r °ai £ ~ °a.2 ™p 



(from i = 0 to |=1) 



,3 



bkr a "- = °b = %o + °bi - °b* 



D 2 y n 

(from £ = 0 to * = 1) 
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4. The wave length of the elastic surface .— We have already- 
stressed .the Importance of the tests on the deformations in in- 
cipient wrinkling, so as to gain at least an approximate con- 
ception of the feasibility of the orthogonal-anisotropic plate 
(with evenly distributed stiffness). . • : 

In this the chief interest centers in the length I of a 
half wave of the wrinkle measured in the direction of the x-axis, 
which must be greater' than the spading of the stiffeners, or' the 
assumption of a mean value for the bending stiffness is no 

longer valid.. The length of the .half-wave of the elastic sur- 
face (measured along the plate length) which begins to form un- 
der critical shear can be computed for. the ideal orthogonal- 
anisotropic plate with one of the two equations: 



l kr '= c a' a / ^ = o a »- a. J p (4a) 



l kr=Cb r a /2L = c b » a y> (4b) 



The absolute factors c a ! a,nd c^' are, like c a and c^, depend- 
ent on # (formula 3) and can be read from Figures 6a and 6b. 
(Factor c a » is used to determine (for # > l) the half-wave 
length forming on the elastic surface at wrinkling (lfc X = o a » 
S'*/p) i Figure 6a; and c D ' "(for & < 1) to determine the half- 
w a ve. length at wrinkling (lfcr-= °b l a<\A0> Figure 6b.) When 
£ > 1', equation (4a), when £ < 1, equation (4b) should be used. 
The approximate formulas for the wave lengths with constants 
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c a'o = 2 '- 05 and 1» 38 > ' Pb'o = 1 - 92 and 1 « 16 

c a \ = 0.44 « 0.28, o- b [ 1 = 0.57 » 0.9 

c-o' 2 =0 . « 0.4 



' D 2 

are 

' n t — o" 1 + o ' 2: " (from i = 0 toi = l) 

c a ~ G a o T u a # v # 5 7 

<V = c b'o + CTd'i ^ - °b'2 ( 11 * = 0 n * = 1). 

for elastically and rigidly constrained support. 

5. Exampl e.- A long, flat, rectangular duralumin plate of 
2a = 100 cm width and 6 = 0.05 cm thickness, like the cross 
section shown in Figure 7 (along x-axis, plate length) is equip- 
ped with riveted angle plates which run parallel to the y-axis 
(plate width direction). The even spacing (distance at center) 
of the stiff eners is d ■ = 4.0 cm. No stiff eners are provided 
in the x-axis. 

We disregard the transverse elongation for the present, so 
that Poisson's ratio becomes v = 0; of course, this does not 
hold good when the shear modulus of the material is G = \ E; 
E = Young's modulus. 

Owing to the difficulty of correct interpretation of the 
elastic edge restraint,, it is no more than a rough calculation, 
so that our simplifying assumption is justifiable. Moreover, 
the cross stress in a reinforced plate is decidedly less effec- 
tive than in a strictly homogeneous plate. 

The angle plates have no appreciable effect on bending 
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stiffness D x so that we can -substitute that of the nonrein- 
foroed plate,, which, ignoring .the transverse elongation, yields 

1 12 

_ 1 x 10 s X "Q :Q5 3 
■ "12 

' = 7 x 1.04 kg cm 
E = 700,000 kg/cm 2 .- " . \ .• 

The stiff eners have likewise no marked effect on the dis- 
tortion stiffness,, so that (omitting v x and Vy) : 

D 3 = 2(GJ) xy 
12 

. = 7 x 1.04 kg cm 

The cross-seotional inertia moment of four angle plates with re- 
spect to the x-axis is 

J = 0.665 cm 4 , 

so that the mean bending stiffness D 3 becomes 

% ~ d 12 

, , = 7 x IQ^x 0.665 + ? x 1>04 
= 7 x 1664- kg cm 

According to (3) 



- 7x7 1.04 X 1664 



7 x 1.04 
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The figures, 2a and 2b, yield the corresponding value': 
c af = 8. 25 (free support), . 

and 

c ae = 15.25 (rigid edge restraint), 
so the shear strength computed with (la) "becomes: 

+ a y 7 x ^ 1.04 x 1664 3 

= 6.1 kg/cm (free support), 

7 x V 1. 04 x 1664 3 . 



and 



t kre = 15.25 X 5 q 2 > ' 

= 11.25 kg/cm (rigid restraint). 

To determine the distance .. l kr of the junction lines on the 
elastic surface, we resort to the- factors given in Figures 6a and 
6b. For £ = 40, we substitute the figures for £ = »: 

c a 'f = 2.05 

and 

c a ' e = 1-38 

where, according to (4a), the half-wave length is defined as 



l krf = 2.05 X 50 j 
= 16.2 cm 



We = 1. 38- x 50 y 
= 10.9 cm. 

In the present case the half-wave length exceeds the center 
distance d of the stiff eners; moreover, the distribution of 
the reinforcements is comparably uniform within length I , so 
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that the calculated shears strength- should he quite accurate. 

In a subsequent section (10, page 49) (Figures 15a and 15b), 
we show that a slight deviation (increase or decrease) in half- 
wave length - in a homogeneous or thogarial-ani so tropic plate - 
presupposes an increase in shear load only slightly lower than . 
that of the shear strength. A certain discrepancy from the fig- 
ure given by formula (4-) always will exist, unless the rein- 
forcements are distributed perfectly even?uy and with vanishing- 
ly small spacing. ■' - • : 

naturally, we must not conclude, herefrom that the ' shear 
strength computed by equation (l)- would always be 'a little too 
low. The local reinforcements affect the position of the junc- 
tion lines (at variance with the approximate calculation) so 
that the reinforced' points on the plate are subjected to lower 
form changes than the -.approximation calls for, and the whole 
work of deflection at wrinkling, and thereby the shear strength, 
is perhaps lower than the calculation warrants. 

In order to define the influence -of the reinforcing plates 
on the shear strength we examine the nonreinfprced 0.5 mm .thick 
duralumin plate (plate strip without angle plate, 2a = 100 cm). 
Again disregarding the .transverse elongation,' the plate stiff- 
ness is^ written : as 

TV = D = D =• § i ■ ■ . - . 

= 7 x 1.04 kg cm 
Including the transverse elongation, we have instead of E 6 3 /l2, 
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for D (Di , D 2 , D 3 ) : 

n = E §3 4 

i2 (i - v*y * 



The D value now is multiplied by a factor 

1 



g v = 



1 - u 3 



Poisson's ratio v for metals is 

t) = 0.25 to 0.3. 
This corresponds to factor 

Cu .= 1.07 to ,1.10. 
With = 1.07. for duralumin, -we. obtain* 
D = 1.07 x- 7 x 1.04 = 7. 8. kg cm. 
Sinoe . $ = 1, we have- 

II o. 

t krf = 15.165 x = 0.041 kg/cm {freely supported) 

50 

t kre = 32.15 x = 0.069 kg/cm (rigidly restrained) 

50 

T krf = §^ ' = 0.82 kg/cm 
T kre = . = 1-38 

Y 

; krf = 2 - 49 x 50 / f 1 ! = 124 * 5 cra 
l kre = 1«66 x 50 = 83.0 cm. 

*3ending and distortion tests on duralumin plates yielded as high 
as v = 0.24. (Compare Z. f. Tech. Physik, Vol. 8, 1937, pp. 355- 
359: M. Bergstrasser , "Determination of the Two Elastic Constants 
in Platelike Bodies.") The equivalent of v = 0.24 would be 
Ov = 1.06. 
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For v = 0, the figures for t^ r and. T kr would be 7$ 
lower; would remain unchanged. 

The critical shear stress' of the nonreinf orced strip lying 
between the angle plates oan be calculated- in the same manner 
(Fig. '7) . • The half -width of these strips- is less than 1 cm. 

Calculating without restraint, we obtain a critical shear 
stress which in any case is higher than 

'.' ty T - 13.165 x 

1.0 s 

• = 103 kg/cm . . 

As is to be expected, this figure is far above that found '-for 
the whole reinforced plate. .. • 

The conformity of very thin isotropic homogeneous plates 
with our assumptions set .up . in the development of the formulas - 
plate thickness small compared to length and width, and small 
elastic deflection relative to plate thickness - may lead, one 
to conclude that the computed elastic surf ace would appear par- 
ticularly well in tests with very thin plates. This, however, 
is not always true, because such plates cannot be manufactured 
absolutely flat and show small bulgee which, owing to the thin- 
ness of the material-^ must 'be taken into consideration. 
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III. Investigation 

6. The differential equation of the elastic surface of an 
or thogonal-ani so tr op ic pla t e . ~ Let 

p = load per unit of plate area acting perpendicularly 
to the median plane; 

My and M x = bending moments per unit of length acting on 
areas of intersection 6 dx and 8 dy; 

H xy = H yx = H = shear moments per unit of length acting on 

the. areas of intersection. 

The equation of equilibrium between load p and moments 

M x , My, and H then reads: 

Nq assumptions as to the elastic properties of the plate 
were made, and the equation applies equally ' to a plate with any 
elastic property or 'with any elongation in plate area. 

Each of these moments, M x , My, and H induces certain elas- 
tic deformations, .but of which only the deformations w (in the 
direction of the z axis, positive, downward) hold any particular 
interest for us. With only a constant bending moment M x (that 
is, My = H '= 0) acting in an arbitrarily long part of the 
plate, the result - when all points of the plate are identical 
in elastic property, i.e., the plate is homogeneous - is 

1) a constant curvature (about y-axis) along x~axis; 

2) a certain (lower) curvature (about x-axis) in the y-direction, 
due to the transverse elongation.- So in order to describe 
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this form change, we naturally mus*t- know the two necessary quan- 
tities on which the elastic properties of the material depend 
, (constant in homogeneous materials). This applies to moment My 
also. With only K (that is, M x = M y = 0) , active, a distor- 
tion results. H acts on all four intersections (5 dx and 
5 dy) of the assumedly cut-out. plate elements, which, to de- 
scribe this attitude of deformation, requires only one elastic 
quantity (a constant in homogeneous material). 

The special (orthogonal) anisotropy of the plate is inferred 
from the absence of distortion by the bending moments M x or My, 
and from the absence of deflection by distortion moment H with 
respect to the x, y axes. Owing to the limitation to small de- 
flections w, the ensuing small quantities are negligible and 

the curvature is expressed as ££-1 and , and the 'distor- 

ts as By 

tion as JLJL . 
• 3x By 

How the described elastic properties of the orthogonal- 
anisotropic plate can be written as 



'x 



= — •■- . ■ 1 ly T —. — «~— — mv 

a x 2 (EJL X <3J) v > 



a 2 w „ .. . i „ M * 



37 - (EJ) y ^ + Tif)I Mx 



g 2 w_ 1 



Bx 3y 3 (GJ) 



H 



xy 



(6) 



with -U x and Uy as certain ideal Poisson's ratios. Resolving 
these equations according to M x , My and H, we have 
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(SJ) X w + ^ 9 3 w 

T ' 9 y 

3 3 w 



H x = 



1 - % Uy ^3 X 



M y =•- 



(EJ), 



1 - v x v y v 9 y 



9 2 w 

— + % 



3 x' 



) 
) 



(7) 



Bx dy 



The five elastic quantities (EJ) X , (EJ) y , (GJ) X y>. u x . and 
u y now are ...defined by formulas . (8) and (7); the elastic quanti- 
ties need. not even be constants; they might also be functions of 
the locus (nonhomogeneous plate). Accordingly, these elastic 
quantities must be defined by experiments in bending and distor- 
tion in ■which. only one constant bending, and distortion moment is 
active. Poisson's ratio p x and p y is generally low enough to 
be set at zero for the conventional plate reinforcements. Be- 
ginning, as in the case . of the flat, homogeneous plate, with 
simple load cases (tension and shear loading) and deriving equa- 
tions (6) by applying some known theorem (Hooke, Bernoulli), the 
five elastic quantities can be determined -from the tension and 
shear load attitudes. But inasmuch as this does not enter into 
the question in the majority of the- systems considered here, we 
assume these elastic quantities defined by (.6) and (7). 

The work of form change of an orthogonal-aniso tropic plate, 
loaded with moments ' M x and M y is unaffected by the sequence 
in which M x and My are used, -and the elastic quantities yield, 
according to Maxwell and Betti, the following relation, 
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Vy (EJ) X = Px (Ej)y 

Thus the orthogonal anisotropy of a thin plate is charac- 
terized "by four elastic quantities, the case of general anisot- 
ropy of a thin plate "by six, and the three-dimensional "body "by 
21- quantities. 

Putting, briefly, 

_ ' (EJ) X 
■ 1 - v v u 



D 2 = 



x "y 

(Ej)y 

1 - v x iy 
(EJ) 



2_ 



2 D 3 = tv 

' 1 ~ U x Vy 



(8) 



and inserting formulas (7) in formula (5) the differential equa- 
tion of the elastic surface of a thin orthogonal-anisotropic 
plate reads 



3 W + p d 
8 x 4 3 c)x 2 Sy s 



Pi 4 w 

9 y 4 



p = 0 



(9) 



7. St ability equation, edge eq uat ions an d theor em for re- 
solving the diff erenti al equ ation.- Subjecting the edges of an 
elastic, orthogonal-anisotropic strip of infinite length and of 
2a width to an evenly distributed shear t, produces at first 
a pure shear stress attitude in the plate, but as soon as t 
reaches a certain magnitude, the plate begins to buckle. This 
limiting case of the stability is given for the lowest load 
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t (t]j r ), and equation (9| is complied with when we set 



-Vi + at if-,- 

. ax By 



and (9) becomes 



Dl 9!l + 2D 3 -^JL + D 2 it? + 2t = o (.10) 

1 3x 4 3 Bx 2 3y 2 2 By 4 3x By 

Since the strip is of infinite length it is necessary that 
the solution be periodic in the longitudinal direction (x) . 
According to' Southwell and Skan, we wrire 

iK-S r=4 iX:jX 
w = e a 2 C r e a . " - (11) 
r=2 

where k is to be real, and the quantities X real or complex. 
This equation, written in (10) then yields 

i k | r= * i^rj _ ' K 4 k b X* 

r=i a a 3 a 2 

X * k X ri 

3 a* a a 

Now we add the real quantities (Compare formula (3)) 

|i = 5a, p = 2x a nd t, = - ~ 

and we arrive' at the characteristic equation with the unknown 
factor X: 

X 4 + 3 u k 2 X s - 3 t x k X + p k 4 = o (13) 
With X 1} X 2j X 3 , and X 4 as the roots of this formula, we have 



(X - X x ) (X - X 2 ) (X - X 3 ) (X - XJ = 0 (13a) 
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as 



Since the term with. A. 3 does not appear in. (12), as a compari- 
son of the coefficients in (12) and (12a) shows, the four roots 
must comply with ■ . ■ , 

X x + Xs + X 3 + A4 = 0, 

and we express them as 

X x = a (+ 1 + n) 
X 2 = a (+ l - n) 
X 3 = a (- 1 + m) 
X 4 = a (- 1 - m) 



(13) 



The coefficients in (12) (+ 3 M- « a , - 3 ^ k' and + P k*.) 
are real; if this equation has a complex instead of real root, 
the conjugate-complex number is a root also. The result is 
that m and n are either real or purely >■ imaginary (that is, in 
no case somplex), and a always real. 

A comparison of the factors in (12) with those in (12a) 



'yields for X, that is,- for a, m and n 



D 



a 



) a 3 (- 3 - n 3 - m 3 ) = 2 U ^ = 2 * 2 



D, 



• b) a* (1 - n s ) (1 - m 3 ) = P ^ = ^ 

^2 



c) cc 3 ( n 3 - m s ') 



= T, K = 



t a 3 
D 2 



(M) 



and further relations between" a, m, and 11 furnish the edge 
equations .„.'._ 

The plate strip is supported on the edges y = + a and 

y = - a, and we assume that there is no deflection on these 
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edges. How, the first two equations are 

a) (w) y=+a = 0^ 

b) (w) y= _ a = 0^ 



(15) 



This conforms to the supposition of bending-resistant edge sup- 
ports.* 

Moreover, we must determine the effect of the type of sup- 
port on the tangent slope — of the elastic surface at the 

edges. If the support is : on knife-edges the plate is n:>t pre- 

♦The inclusion of bending of the longitudinal support could be 
accomplished by means of one of the following equations for the 
elastic line of the edge support with (E'J^) bending stiffness, 
namely ■ _ 

(s ' J *> (MO + - * 

°"- y y= ± a 

In the first formula the bending moment at point x of the edge 
support is calculated as the sum of all turning moments of the 
cross stress Q y and of the turning moment H of the plate at 
the cut-off part; that is, when the girder is nowhere supported: 

Mb = / [Q Y (x - I) + H] d. t 

In the second equation the loading p x per unit of length, com- 
prises the cross stress Q, y and the derivative of H, that is: 

The edge girder need no support because, due to the periodic 
character of the solution, the edge stresses are in equilibrium. 
If, however, a support is provided, the latter must be evenly 
spaced (l) so as to conform to the periodic solution (equation 
(ll) ) and, in addition, the solution is generally narrowed be- 
cause of the added formula, 



< w )y=±a, x=nl = 0 
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vented from such a slope Of -its elastic surface. But as a rule, 
the plate is so fastened to the edge support that, "by perfect 
attachment, this support distorts "to an amount equivalent to 
angle -~ of the elastic" surface' on the bearing. If the edge 
support offers an elastic" resistance against this distortion, 
it produces bending moments My in the plate and distortion mo- 
ments (per unit, length) 'in the edge support of the 'same size as 
the plate bending moments, and the restraint is elastic. How 
it may be inferred that the increase in distortion angle (per 
unit length of support axis-,- or ; .,of the edge line) of the edge 
support is at every point proportional to the (inside) turning 
moment of the edge support. This corresponds, strictly speaking, 
to a support of constant,,' round ..cross, section with uniform dis- 
tribution of the turning moments over the cross sectional per- 
iphery. But in first approximation it could equally be applied 
to other than round cross sections. 

At the edge y = + a the' distortion angle at any point is 

. . ... Vd y/ y^a ■ 

and the increase in. distortion. angle per unit, length is 

d 6 = / d 2 r; \ ■ .. 
d x V8x dy/ 

y=+a 

the distortion load; m^ applied: to the support per unit length 
now equals the bending moment (My)y=+ a of the plate; hence, 
according to (7) and (8): 
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The inside .distortion moment acting at point x is found by in- 
tegrating over the .distortion load m&, assuming the edge sup- 
port, to be unprotected against distortion. 
We write: 

c MIf) ]dx) * 

vay / y=+a v « / y ^ +a 

With (G'Jcl) as distortion resistance, we have 

d5 = M d . 
dx W Jd)' ' 

and, by inserting the corresponding values 

x 

/ d S W' \ - 1 . f fr, ( f w \ . , 

V9x oy/ y=+a ~ (G' J d ) lD * \ a y 3 / y=+a 



VOX -y=+ a 



For edge y = - a the conformal equation (since distortion 

3y 



-.3 

' l£LaB a corresponding distortion moment m^ = - (My)y = _ a 



is 



*The integration can be avoided, as in the preceding footnote, 
by using another (equivalent) form of torsion balance of the 
edge support, namely, 



instead as above x 

(G» Jd) f§ = % = - / md dx. 
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(_£jL> = + 1_ / [D ' + 

Vas (Gi J * } ■ SyS/ y=-a " 

Now we insert (ll) in (15a) to (I5d) and have.*. 

. x _« 

(a,b) e 1Ka r E 0* e ±iX * = 0 • 
r=i 

(o,d) \ e 1 ^ X T C r X r e ±iXr 
a^ r=i 



x 

(&'' J d ) i k a ~ u ^ 2 r=i' 



= ± —1- - e a [D 2 S b r - V e* 1 ^ + 



+ % D s ^ *|* 0 r e ±iX r] . 

In the terms with double signs,, the' upper sign (equations 
a and c) is valid once only in each equation, and then the lower 
sign for equations b and.d. The last term in. (c,d) disappears 
according to equation (a,b). The relation for % 8 



e 1 r = cos \p ± i sin Xp 



and we abbreviate 



D 2 a 



(G» Jd) K * 

Now the equations (a,b) and (c,d) become: 
r=* r=* 



(a,b) .... 0 r cos Xj, ± i 0 r sin Xj, = 0 

/ r=* . r=*. ' ^ . 

(c,d) ....( 2 0 r X r cos X,, + e 2 0 r X r s sin \ i 
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Setting the real and the imaginary parts of this equation to 
zero, the four-edge equations now read 



r= 4 



a 



) Or c ° s X r 



b) X G r sin. X r 



= 0 
= 0 



r=i 
r=4 



(16) 



°) r ^ x °r \r cos ^r + e r 5 °r V* sin ^r = 0 

r= * r= 
d) 2 C r -Xj, sin X r - c 2 0 r X^, 2 cos X r = 0 



These are four homogeneous equations with the unknown 0 
factors. To assure ' nondisappearing 0 values, the determinant 
of the denominator must disappear. This yields a transcendent 
equation with the four X roots; by introducing (13) we estab- 
lish a relation between a, n, m, and e. 

e depends on the form of the edge supports and may, accord- 
ing to (G' J&) assume any value between zero and infinity. 
To simplify the calculation, we use only these two limiting val- 
ues of e in the subsequent examination. 

For 

€ = 0, 

it corresponds to an infinitely high (&' J^) , value, that is, 
an edge rigid in distortion; this is the case of (rigid restraint 
(case "e"), where .the angle of tangent slope |2 of the elastic 
surface equals zero at the edge. 

Now we divide (16c) and (l6d) by e, and make 



i = 0, that is, e = <*>, 



(17) 
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which corresponds to the knife-edge support (free support, case 
"f") of the plate edge, and where bending moment (My)y = ± a dis- 
appears. 

When we insert the given values, and determinant of the de- 
nominator becomes zero, the relations between a, n,, and. m, for 
cases "f" and "e" are: 

, > 
f) ' ; 8 n'm [cos 2 n a cos 2 m a - cos 4 aj 

-'•[4 (n 2 + m 2 ) - (n 2 - m 2 ) 3 ] sin 2 n a sin 2 m a = 0 

(by free support) 
e) 2 n m C cos 2 n a cos 2 ma- cos 4 a] 
- [4 - n 2 - m 2 ] sin 2 n a sin 2 m a = 0 

. (by rigid restraint) 

8. The critical shear stres s.- The equation (14) combined 
with one of (17) yield the relations between D x , D 2 , D 3 , a 
and that shear stress t, at which the pure shear stress ceases 
to present a stable attitude of equilibrium, but where the equi- 
librium equations yield a possible state of deflection conforma- 
bly to equation (ll). Hereby it is assumed that the stresses 
consistent with this attitude do not exceed the proportional 
elastic limit.. 

The. four equations - (14a, b, c), (I7f) and (l7e), respec- 
tively, contain five unknown factors: k, a, n, m, and t. The 
k value determines, in accordance with equation (ll), the wave 
length of the elastic surface in the direction of axis x. 

In a definite system with definite Di , D 2 , D 3 , and a, we 
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may at first assume an arbitrary wave length., or a suitable k , 
and then determine the respective t from the four equations. 
The desired shear strength t^- r is the lowest t value at 
which the plate commences to wrinkle, or in other words, where 
a function w (equation 11) different from zero, prevails. 

To extend the applicability of the systems considered here 
(for arbitrary , D 2 , D 3 , and a) we transform the equations. - 
We introduce in (14) the factor £ for. the orthogonal-aniaotrop- 
ic plate defined from equation (3) 

D 3 M- 

which may assume any real value from aero to infinity. Then 
we introduce the parameters 

a) 
b) 

Now we divide- (17) into a transcendent function cp equivalent 
for ca.se s f and e, which contains a, m, and n; and into 
an algebraic function t different for cases f and e, and 
which contains n and m • only. 
The result is 



= Ja 
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a) 
b) 
c) 



_ - ta and. « £ b 

(l ~ n 3 ) (1 - m?) = t*- and = * a 



JbI D 2 3 a 3 (n 3 - m 3 ) 



a^ 



' D 2 D 3 as (n 3 - m s ) 
and = - ■ a -_- 3 x 



.2- 



(19) 



With Bergmann and Reissner designations (19c) becomes: 



t a 3 _ a 3 (n 2 - , m 3 

K 



1 - Vn a 3 (n 3 ~ m 3 ). 

and ~/jT (nL" 

: ■ V % . 

f) cp = tyf (by free support), 1 
e) 9 = ' ty e (by rigid restraint); 



it denotes 



cp = 



sin 3 n a sin 3 m a 

cos 2 n a cos 2 m a - cos 4 ct 

= w . tan 2 n a tan 2 m a 



if we make 



co =• 



1 - 



cos 4 a 



and 



cos 2 n ex cos 2 m ct 
8 n m 



^ - = , 

1 ■ 4' (n 3 + m 3 ) - (n 3 - m 3 ) 3 

^ • • _ ■■' 2 n m 

e . -.4 - (n 3 + m 3 ) 



(20) 
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This is simplified with 



a ) n - aS (n s - m 3 ) 

a) ° a 7Z 



= a3 ~ m? ) 
= CoV^ (See equation 18b), 



(21) 



and the shear strength t of equation (19c) now yields 



a) t = 0 a 



b) t = o b 



a" 



7d s D, 



a' 



" 3 



For the case that t represents the lowest value actually 
obtainable, i.e., shear strength t kr , we have 



/ D D 3 
a) ■ t Vr = c ; ^ 1 a 



a 



a 



c kr 



= c 



a 



b) t kr = c 



J kr 



a" 



= ob. 



when 



and 



a) c a = [ a3 (n! =_slj] 

b) c ^ [ aS (nlr )] 

v St 



mm 



mm 



(21') 



These are values which, as we shall show, may be represented as 
simple functions of & . 

Ngw we proceed as follows: We determine n and m as terms 
of and of £ a or £b » conformally to (19a) and (19b); 
then we compute the corresponding a from (20f) or (30e), and 
0 a and Ob from (21). 0 a and % assume, as £ a and ^ func- 



N . A . C . A . Technical Memorandum No . ■ 602 



35 



tions, certain minimums for each & (value). These minimums 
represent the desired coefficients c a and 0-5. . 

Now we ascertain whether C a and 0^ and thereby t as 
functions of & and ^ (and £5) .are so. formed that these func- 
tions represent a solution of equation (ll) when $ assumes any 
value between zero and infinity; then we must find the range of 
£ a and £b as functions 'and their minimum value. 

We begin with the partial "functions which comprise function 
0 a (and 0^) - equations (21a) and (21b). 

From : (l9a) and (19b) we. calculate n and m as well as 
(if + m 2 ),- (n 2 - m 2 ), and n mi . . 

C, 



n 2 + m 2 



and 



= - (1 + Cb) 



_ + 



2 + 



- c 



a 



and 



= + /(2 + £ b f _ £ 2 S b 2 



n 



n s + m 3 + nj 



m" 



= ± 



and 



y- (1 + u +7 (2 + ^b) 3 - * s £*> s 



m = ± 



= ± 



n 2 + m 3 n 2 - m 3 - 



. 2 



1 + 



£a\ _ /(?+ laX 5 ~ t 2 



(22) 
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and 



n m 



/ ng + nP^ a / a? - nPN a 
V V 2 V ~ V 2 / 



a 3 Ca B 



- y £a-.. 5 



and 



= J £ 3 £b 3 ~ 2 £b - 3 



According to assumption (equations 11 and 13) k and a, 
as well as # (equation 13) are real; thus it follows from (18) 
that £ a and ^ are always, real and positive. But, according 
to (22), (n 3 + m 2 ) are always real (and negative). So n and 
m are either real or purely imaginary, but. never complex (Com- 
pare equation 13); this is feasible only when (n 2 - m s ) is 
real also; or, in other words, (©ompare equation -22) when 

t 2 



( v 2 T) a ' 



Then we can write 



3 + f = i £ 



For * = 1 this equation is complied with for every posi- 
tive j^, while & > 1 yields 



and 



f ^ 2 g / _ 2 

^ (for 

^ _ a - 1 



>D 



(23) 



For these £ a and - £b values, 
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Figures 8a and 8"o show n and m as ordinates plotted 
against the abscissas £ a and C.b with S- as parameter. In it- 
self (aside from the limiting cases) one of these two would suf- 
fice ( £ a or £> D as abscissa). But for & > 1 the numerical 
calculation is more simple when £ a is used as variable, but 
for # < 1, is preferable. Accordingly, the curves are 

more comprehensive when we select two different variables £ a 
and .'k for different For & =::1, ■ £ a = ' and the shape 

of the curves' is the same (Compare equations .18 'and 22). 

Within the range of £ a and ^ considered here, '"and whose 
upper limit for & > 1 is -given in equation (23), m is always 
imaginary (equation 22). For imaginary values of m and n we 
write, for short: " .. : .;' ' /./../•. 



m = B N 



- " I 



(24) 



For t a = 0 (and = 6)' m has (independent of &) 
the value m = y~3~. 

When £ a ( ) increases,, m increases at first (except 
for £ = oo ) . For 1 ^ & = », m ..becomes by 

f - 4 ^ 



maximum 
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and then drops to a minimum: 

m = 

which is reached by the limit value £ a in (33). For 
0=^1, m tends uniformly to infinity. ' For £ a = 0 ({■£, = 0) 
n is maximum 

n max = + 1 (independent of £ ) . 
As £ a (£b) increases, n remains at first real and then 
drops evenly to a minimum n = 0, which with a perpendicular 
curve tangent is reached at 

r i + 7 l + 3 £ 3 
^a j 

and 



1 + J 1 + 3 & 



(25) 



When £ a (£b) exceeds this value, jj_ "becomes imaginary. For the 
highest £ a value considered here (equation 23) curve n with 
perpendicular tangent changes into curve m. Likewise, of im- 
portance are functions 



and 



a) X^^St 



•D 



(26) 



illustrated in Figures (9a) and (9b). Both functions have an 
infinite value when £ a = 0 and ^ = 0. For £ > 1, / a 
turns uniformly to sero, which it reaches with perpendicular tan- 
gent by -• 



N.A.C.A. Teclmical Memorandum Uo. 602 



39 



£a ~ t (Compare equation '23). 

For # < 1, %b reaches a minimum 

• ^bmin = + J 1 - * a 

at •'• . • 

£b - ~ 1 :• (260 

y 1 - # 2 

Beginning at this minimum goes evenly to infinity when 

^ increases. & = 1 yields a. horizontal asymptote with ordi- 
nate 4. So.,., while X a ana ^ "b decrease when C a ( £- 0 .) increases 
beginning at £ a = 0 and (£t> = 0), a. and a 3 increase .and ., 
approach infinity, so that the product as % & (and as t %.^ 
(Compare equation 21) , assumes a minimum value. 

Every d- value and. every £ a (£t,) value within the range 
defined in (23) has a certain pair of n and m .values (See 
equation 8) and, if beginning at £ a = 0 ( $b = 0) , n is at 
first real and m imaginary. This is why the letter i occurs 
in the numerator of <p and ^ in (20), and which is expressed by 

<P _ sin 2 n a sin 2 m a 

i GQ~2~n"lx~ c^~2~^T~cT~'5'5i"~4 d 

= tan 2 n a tan 2 m a 00 

The function cp/i* is shown plotted in Figure 10 for an arbi- 

trarily cho sen pair (according to equation 22) 11 (= 0.87) ' and 

*The aim-is to express that m as imaginary. If ~m and h "are~~ 
imaginary the' expression cp/is (instead of ~<P) is used. \j/ /i 
and t/i 2 are used in the same • manner . 
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m (= 1.73) which corresponds ' to ■& = 10; £ a = 0.98; the figure 
also includes tan 2 n a and tf/i and i e /i plotted against 
a (See equation 30). In a subsequent calculation of the criti- 
cal shear stress \^e intend to show that n is never much lower 
than "one"; and even for materially lower n the course of 
function <p/i is not much different from that on Figure 10. 
The zero points of the function are usually at a values which 
form an integral multiple of Tr/3n, while the functions at 
those a values inhere the denominator disappears, have perpen- 
dicular asymptotes and change from + 00 to ~ 00 . With a = o 
denominator and numerator disappear; double differentiation of 
numerator and denominator of function <P yields 

/ \ _ 3 n m 

^ ; a=o 4 - (n 2 + m 2 ) 

This, according to (20) is the value for t e , so that 



comes evident for the reason that 

Cp(Oc) = 9 (-a) 

so that the function is symmetrical with respect to the ordinate 
axis. From a = 0, 9 goes evenly toward infinity. When 
1 > n > 2/3 there is an asymptote for a value a, which is 

1 TT 

somewhat higher than ~ — ; in any case there is an asymptote 
near this value and another asymptote near ~ — , and between 




By a = 0 function <p/i has a horizontal tangent as be- 



4 n 
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the two <P intersects the axis of the abscissa when a = i _ . 



Each cp/i function with a pair of n, m values has certain 
tyf/i and te/i values (equation 30.) which are shown in Figures 
12a and 12b plotted against ■& and £a(£b)> tf/i is always neg- 
ative, .and t e /i positive. In Figure 10, tf/i and t e /i, 
being independent from a, are parallels of the axis of abscissa, 
and at a = 0 .the straight t. e /i, is a tangent of function cp . . 
The intersections of the straight with curve . <P yield a as 
solution of equation (20). .For £a = 0 ( £h = 0), tf/i and ^ e /i 
yield. an absolute maximum, which is equivalent for every & 
value. 

For other possible pairs of .n.and m then in Figure 10, 
the same conditions prevs.il, when a is used as multiple of TT/n; 
the position of the perpendicular asymptote and the ordinate of 
(q> ) a=0 is slightly shifted. For n values near zero the 
straights ' tyf/i and %/i are very close to the axis of the 
abscissa, and coincide, like cp, with this axis -when n = 0. 
Aside from this special case the solutions of equations cp/i = 
* tf/i yield a -values which are slightly-lower than the integral 
multiple of V&hy while qp/i = te/i yields, aside from a = 0, 
other a- values which are slightly higher than the, integral mul- 
tiple of ■ Tr/2n. .■" 

Being able to compute oc by approximation,., facilitates the 
solution. According to (20)' <p/l can be decomposed into 



%' = tan 



2 n' a tan 2 m a ax 



0) = 



1 



1 



cos 4 a 



cos 2 n a cos 2 m a 
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tan 2 m a lies between 0 and 1; when 2 m a increases it ap- 
proaches "one" very quickly. Since m here is greater than n 

(in most cases twice as great as n) , tan 2 m a is not much 

tan 

smaller than "one" for a * tt/3 n (for example,/ 2 tt = 0.999993). 
The <d value for a = l/8 tt, 3/8 tt, 5/8 tt, etc. is "one" (un- 
less n = 1), because then cos 4a = 0. By a = 1/4 x tt /n, 
3/4 x Tr/n, 5/4 x rr/n, etc., to = 0, because cos 2 n oc = 0. 
For a values which satisfy equation cos 4- a = eoc 2 na cos 
2 m a j, co = ±», because the denominator of w disappears. 
On these points we find the asymptotes of curve cp/i. Figure 11 
shows w plotted against the cp/i curve of Figure 10. Near 
a = tt/3 n, we approximate cp/i as 

%' « tan 3 n a . 

Having defined an approximate .value' for a we. compute 
go tan 2 m a and estimate the error of the approximation. The 
dotted' line in Figure 10 represents the approximation curve 
tan 2 n ex. 

. Of all feasible solutions of equation cp =• i|r , the solutions 
a near tt/2 n represent only those values which yield the 
critical shear stress. Presuming this assertion correct, the 
functions 0 a = a 2 X a and 0^ = a 3 become minimum. 

When : f *' > 1, • then a and as are finite for £ a ' = 0, and 
increase evenly when £ a increases ; for n becomes smaller as 
L increases (Compare Figure 8a). Since the changes of t f 
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and. i|/ e in equation cp/i = \|//i (when £ a changes) have a rel- 
atively small effect on the a solutions, the a values are ap- 
proximately inversely proportional to n and thus increase 
with £ a ; when n goes to zero, a becomes- infinite. But Xa 
decreases evenly in this same range from Xa = 00 ( for ^a = 0 
and assumes for £ a , where n = 0 (See equation 25) a finite, 
real value. For £ a = 0 and for £ a conformal to (35), the 
product of a. 2 and X a becomes also infinite, while between it 
it assumes finite values. 

For < 1 the Xb function becomes infinite at ^ = 0 
■and ' £b = 03 ; between these two is the minimum indicated by (26*) 
a becomes finite for £- D = 0. Aside from the limiting case 
& = 0, n decreases as £b increases; when. S-j-, approaches 
the value given by (35) - where n = 0 - then a becomes infi- 
nite. Between, this value of ^ . and £-5 = 0 the product of 
functions a 3 and Xb ' a minimum value. 

■ That' there is a minimum for a 2 X"b even if & = 0, was 
proved by Bergmann and Reissner. 

When n is very small, a approaches infinity. When 
n ■= 0, . any 0- satisfies equation <p/i - t/i. 

When £ a (£b) increases beyond the value given in (25), for 
which : ".-n = 0, then n becomes imaginary. Since in this range 
(n 3 m 2 ) and Xa and Xb (26) decreas'e still further, it 
seems feasible that the minimum of 0a(°b) (equation 31) lies in 
just this range. With (24), equation (20) now becomes 
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and =h 

i 3 IT i 3 

cp 

■p- = to t an 2 n a t a n 2 m a 

(o = i_JL_ 

1 c os 4 a 

cos 2 n a cos 2 m a 

Jj£ = 8 n m 

i 2 ~ 4 (n 3 + m 2 ) + (m 2 - n 2 f 

t P 2 n m 

i 2 4 + (n 2 + m 3 ) 

Now to and the product tan 2 n a tan 2 m a are always 
positive, consequently cp/i 3 is positive also. But ^f/i 2 is 
always negative. In the case of free support this range of 
^a^b) yields no solution for a, except when n = 0 (n - 0). 
Otherwise a solution would he practicable only in the case of 
rigid restraint, "because ^f e /i 3 is always positive. 

The shape of curve cp/i 3 with respect to a is like that 
of curve cp/i; Figure 13 shows g?/i 3 and tan 2 n a tan 2 m a 
and to plotted against a for n = 0.62, m = 1.58 (equivalent 
to & = 5, £ a = 2.2)) and for n = 1; m = 1 (equivalent to 
# = «, Ca = 2.0), The straight line te/i S is again (as te/i 
by cp/i) the tangent to curve cp/i 2 at point a = 0. While 
cp/i differs only slightly from tan 2 n a when a increases, 
the cp/i 2 function proceeds like the function tan 2 n a tan 2 
ma (in Figure 13) , that is, approaches "one" asymptotically. 
The difference between cp/i 2 and function tan 2 n a tan 3 m a 
is expressed in to, which is also shown in Figure 13. When n 
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becomes very small, <p/i 2 again (like <P/i) approaches the 
axis of the abscissa, and coincides with it when n = 0 (n = 0). 
By the maximum value of £ a for n and m as given by (23), 
we obtain value "one" when £ = <©. All cp/x 2 curves for 
d- = oo lie between this upper limit curve and the axis of ab- 
scissa. For smaller # the limit curve (where n = m = 1) 
shifts upward, n and m become infinite when # decreases and 
C a increases. The limit curve <p/i 2 (where n = m) approaches 
when & decreases toward "one," the straight line with the ordi- 
nate "one." It will be seen that in this range (for any n and 
m) the only solution of the transcendent equation is the value 
a = 0 for every # value. For n = 0 any value a fulfills 
the equation cp/i 2 = t/i 2 (like cp/i = y/i). 

The next problem is to find which of the various solutions 
a yields the minimum c a = cc 3 \ a and c^ = co 2 ^ (equation 21). 
We begin with a = o and case n = 0, where any cc represents 
a solution of <P = t. 

If a = o, we see from (13) that .all four roots N disap- 
pear, i.e., they are of the same magnitude (^ x = \s = ^3 = ^4). 
In that case (ll) does not yield the general integral for re- 
solving the theorem in (ll); instead of the four arbitrary con- 
stants we have only one. With four equal roots X r = \ x the 
solution (conformal to ll) of the differential equation reads: 
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The case of n = 0 is treated in the same manner; it yields 
^•i = ^2, and in place of equation (ll) for w we have 

IkS [ iX.Z iX 3 Z i\ 4 Z ] 

w = e a |e a [ d + G a ~j + e a 0 3 + e a 0 4 

Both cases (a = o) and (n = 0) yield a new theorem for 
function w, and equation; (16) supply new equations whose deter- 
minant of the denominator, when made zero, leads to a relation 
between n, m, and a as expressed in theorem 1 (equation ll) 
by equations (17). The final result is that there is no func- 
tion w - not zero - .in the form called for by two or more 
equal roots X r . (Southwell and Skan made this examination for 
£ = 2, and Bergmann and Reissner for # = 0, in Zeitschrift 
fur Flugtechnik und M§torluf tschiff ahrt, Vol. 20, 1929, pp. 477-8.) 

This leaves us the solutions a which are near the integral 
multiples of Tt/sn and which resuit from an equation cp/i = ty/i 
(where n is real) • 

Since we seek the minimum value of shear* t, that is, of 
0 a and 0-5 (See equation 21), the solution & near Tr/2n is the 
only one to be considered for every S- value. For the case (f) 
of free support a is slightly lower; for (e) rigid restraint 
somewhat higher than tt /2n. 

N©w we can develop a formula, for c a and c-^ (equation 21 r )- 
Within the range considered 

1 

a = H arc tan & 
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can be vexj accurately approximated by introducing for co a 
corresponding approximative function or simply a constant 
(0) « i). 

Writing the approximated a in (21) it is probable that a 
•differentiation yields an analytical function for the minimum 
Ca(°b) with 'variable But this function which the equation 

would (approximately) represent (See Figures 2a and 2b) would be 
very complicated and unsuited to practical use. For that reason 
we preferred to calculate a series of points on this curve and 
to express the result graphically in Figures 2a and 2b. 

9. Calculation of c a and c- D for various £ values.- For 
# > 1 we have a range of £ a limited by £ a = 0 and £ a as 
defined in' equation (25). For $ < 1 the range is between 

= 0' and S *b as defined in (26' ). For any £ a and t-^ values 
within this' range, we have 

^ % - } ft 3 - ~ , n, m and n m (according to eq. 22) 

X ( » » » 26) 

t and a . ( '» » » 20) 

0 a « Ot ( » " » 21) 

The sign of n and m is, according to (22), either posi- 
tive or negative, but we confine ourselves to the positive, be- 
cause the negative sign yields the same result. For, changing 
the sign by m(m) or n that for cp(<p/i) and ^("^/i) changes 
also, so that (20), (21) and (26) are not altered. The sign for 
a , which likewise can be positive or negative (Compare (20) ) 
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lias no effect on the final result, because it is affected only 
by the direction of t (Compare (14c), which may be arbitrarily 
assumed. 

The solution of (20) is approximated with 0) = 1 [that is, 
a = (-£• n)'arc tan the solution a being inserted only near 

Tr/2'n. This is a very satisfactory approximation for the case of 
rigid restraint. 

By £ a = 0" (t^ = 0) we find n - 1 and m = Jz~ for' all 
& values; and as a minimum we find: 

a = 1.306 (approximated = 1.309) by free support. 

cc'= 1.8330 ( • " - = 1.8323) by rigid restraint. 

and , 

As example, Figure 14 shows the functions t Xa» a > 9> m /i>/ 1/10 
C a plotted against £ a for £ = 5 by free edge support. The 
range of £ a has been enlarged beyond practical use. For 
* = 1 the minimum c a (c^) of C a (0b) is at £ a > 0.7 
(Ck > 0.7), f or £ = 0 at £ b < 1.35 and for S- =. »" a t 
5 a « 1.0. .. .. 

The range of £ a and for practical purposes is: 

^ = 0.7 to 1.0 for # > 1, free support 

C b = 0.7 11 1.35 » * < 1, » « 

C a = 0,9 "l.i « ' * > 1, rigid restraint 

£ b = 0.9 » 2.25 " * < 1, 11 '« 

For * = 0, 1/5, 1/2, 1 and * = 2, 3, 5, 10, 20, 40 and », 
the calculations, are shown in Table I and tabulated in Table II. 
The figures for c a (for £ S l) and c^ (for & ^ l) 
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were rounded off to five units, the calculation was made with 
corresponding accuracy (i.e., the .error for c a and c^ does not 

exceed 0.3 per cent). In Figures 2a and 2b, c a and 0-5 are 
plotted against £ . 

10. Equation and calculation of c a * a nd c p 1 derived for 
the wave length of- the elastic surfac e.- Equation (ll) expresses 
the elastic surface of the strip when wrinkling under shear. 

It contains, aside from a, . the five constants' k } o x , C 3 , 0 3 , 

0' ' 

and 0 4 , of which the four /constants can be computed from the 
linear homogeneous equations (16).. The k constant, defines the 
wave length of the elastic surface in the x-direction. Accord- 
ing to (18) we have 

v 

a a 

There are five constants for every pair of & , £ a and 
# , or, in other words, every pair has a corresponding elas- 

tic surface, which in a homogeneous or thogonal-ani so tropic plate 
would become that of the shear t ^ r . •' But since the usual stiff- 
ened plate never exactly agrees with our assumed ideal system, 
the elastic surface of a p.late when wrinkling deviates more or 
less from the calculation here. 

Since 



k = a JT a /% = i^fi and = a /K = l£k 

* 1 t/P~ b 7 Da f\T 



e a - cos (k + i sin (k 2^; 

\ a/ v a/ 



we infer from equation (ll) that the half-wave I of the elastic 
surface equals half the period of the sine function, that is, 
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So, when we introduce k conformal to (18) the wave length 
is computed as 

4, 



/ 

a) I = 0 a » a J ^ 1 = o a « J^P~ 



b) l= Va y| ™ = Oh' v^" 

where 



(27) 



<V = ^— ; . OV = - C» (Compare eq. 18b) 

If the respective shear equals the shear strength t^, 
we substitute and c a f and c-- D ! for I, C a « and 0- D ' , so 

that (4a) and (4b) yield 



l kr = °a' a ./ ^ = °a' 



The calculation of 0 a l and 013' made by equation (27) for 
various £ a and ' ^ is shown in Table I (columns 9 and 13). 
The values for c a * and c^ ' >are accurately computed to 1 to 3 
per cent. 

Figure 15a shows 0 a m t aP/^/lk Efe 3 \ plotted against 

° a ' = a tfW ' Wlien 1 ^° aI ' chan S es near the minimum (c a ) of 
0 a , the change in t (0 a ) is comparatively slight, particularly 
by free support. 

Figure 15b exhibits % = t ^/J% plotted against 

p ^ 1 

°b = - J\T ■' Tilis ii -'S aiie likewise shows the slight change in 
t (Ofe 1 .-) when I changes near the minimum c^) of 0^. 
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TABLE I 



Calculation of coefficients C a (C^) and C a t (C^t). The figures for the shear strength t kr , i.e. for 

minimum C a and C-5 are tinder scored. 



1 


2 


3 


4 i 5 


6 


7 


8 


9 


10 11 


12 13 




h 


n 


m 

m 5? r* 
— 1 


h 


Free 


Support 


Rigid Restraint 


i 


a 


C -b 


Obi 


te 
~ — 
x 


a 




— 0 


Ac< 
to < 


wording 
aquation 


22 


22 


26 


20 


20 


21 


27 


20 


20 


21 


• 


0 


0.3 
0.5 

1*0 
.1*30 
1 » 32 

1 e«s4 . 

2,00 
2.22 
2.24 




t 

1 


1.898 

2 C 000 
2.236 
2c 3664 
2,3749 


8.400 

7.075 

6.000 

5.7885 

5.7794 


0.4803 
0.4320 
0.3440 
0.3055 
0.3033 


1.3469 
1.3670 
1,4047 
1.4222 

1.4233 


15„24 
13.220 
11.839 
11.7085 
11c 7079 


4.265 

3 «Stj5 

2.238 

1,9374 

1.9212 

1.9053 

1.529 

1.353 

1.047 

0.724 
0.452 
0.367 
0.284 


0.576 

0.5715 

0.559 

0.5292 

0.5225 
0-5219 


1.8325 
1.8304 
1.8256 

1-8144 
1.8116 
1.8113 


28.20 
23.72 
20.04 

18.6227 
18.5894 
18.5889 


3.13 g 
2.43 

1.72 h- 

0 

K 

1.224 | 
1.1639 0 

1.1588 pi 


2,3833 
2.64575 
2.7276 
2.7350 


5.7796 
5.6568*) 
5 . 6646 
5.5659 


0.3010 
0.2405 

0.2078 

0.1443 
0.0739 
0.0249 
0.0143 
0.0072 


1,4244 
1.4528 

1.4684 

1.4990 
1.5340 
1.5583 
1.5636 

| lc566 

1 


11.7081 
11.939 

12.260 

13 0 450 

16.640 

23.88 

28.60 

36,070 


2.26 

2.5 
3.0 
4.0 
8.0 
20.0 
30.0 
50.0 


2.7423 

2.827 
3 t 000 
3.315 
4.359 
6.554 
7.933 
10.149 


5.6674 
5.690 
5.775 
6.000 
7.071 
9.850 
11.690 
14.708 


0.5213 

0.5142 
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0.474 

0*396 

0c285 

0.240 


1.8111 

1.8083 

1.8025 

1.7921 

1.760 

1.7100 

1.6885 


18.5896 

18.62 

18.75 

19.27 

21.89 

28.83 

33.28 

i. - 


1.1539 g, 
1.10 g 
1.01 

0.88 § 
0.63 

0.41 g 
0.34 w 
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TABLE I (cont.) 
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3 


4 


5 


6 


7 


8 1 9. 
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¥e 
i 


a 




j 
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22 

n 


22 


26 


20 


20 


21 


27 


20 
i 


20 


21 


27 


1/5 


1.24 

1.27 

1.28 

1.30 ' 

2.00 

2.04 


0.995 

0.995 

0.995 

0.995 

0.9899 

0.9896 


2.339 
2.352 
2.356 
2.365 
2.642 
2.657 


5.8022 
5.7857 
5.7806 
5.7707 
5.6285 
5.6282*) 


0.3121 
0.3086 
0.3074 
0.3051 
0.2395 


1.4260 
1.4279 
1.4285 
1.4298 
1.4681 


11.798 
11.797 
11.796 
11.797 
12.13 


1.978 
1.953 
1.944 
1.927 
1.513 
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1.833 

1.8302 

1.8302 


19.40 

18.8527 

18.8523 


1.506 
1.214 
1.202 


2.05 
.2,08 


0.9895 
0.9893 


2.6607 
2.6718 


5.6282 
5*6285 


C.521 
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1.0302 
1.8303 


18.8526 
18.8538 


1.199 
1.195 


1/2 


0.4 
1.08 

aao 

1.12 

1.40 

1.50 

1.60 

2.00 

2. SOS 

4.00 

6.00 

8.00 


0.995 

0.9758 

0.9755 

0.9743 

0.9629 

0.9535 

0.9538 

0.936 
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0.6451 
0.4064 


1.947 

2.2610 

2.2697 

2.2780 

2.3932 

2.4329 

2.4718 

2.62 

3.2645 
3.7942 
4.2621 


7.564 

5.836 

5.819 

5.800 

5.624 

5.583 

5.550 

5.475 

5.464*) 

5.653 

6.055 . 

6.477 


0.4564 
0.3304 
0.3277 
0.3249 

0.2335 

0.1260 
0.0710 
0.0340 


1.363 
1.4459 
1.4476 
1.4506 

1.556 

1.860 

2.3800 

3.824 


14.05 
12.199 
12.195 
12.206 

13.24 

19.567 
34.2S8 
94.8 


3.65 
2.091 
2.069 
2.144 

1.427 

0.845 
0.538 
0.448 


0.573 

0.541 

0.5237 
0.5182 
0.5125 
0.4904 

0.408 
0.269 


1.840 

1.864 

1.8818 
1.8881 
1.8951 
1.921 

2.013 
2.658 


25.64 

20.22 

19.915 
IS .903 
19.931 
20.15 

22.85 
42.85 


2,702 

1.605 

1.411 
1.359 
1.311 
1.156 

0.780 
0.482 
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0.73 
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0.8595 
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I 2.1817 
i 2.2044 
! 2.3040 
J 2.3967 
I 2.4211 
| 2.4721 
| 2.5515 
| 2.5809 
2.6832 
2.7404 



10.194 



6.2336 
6.1578 
6.1337 
6.1101 
6.0000 
5.8119 
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23.77 

42.56 
69.20 



2.000 

1.335 

1.000 

0.665 
0.500 



i 0.5736 ; 1.8382 

1.858 
1.903 



2.5829 
2.5147 
2.4928 
2.4713 
2.3690 



0.5122 

0.5098 
0.5086 
0.5074 
0.5049 



34.444 

25.78 

22.75 



1.9527 
1.9580 

1. 9607 
1.9634 
1.9690 ! 

2.079 
2.208 

2.355 

2.694 



22.1600 
22.1528 

22 ,.15 13 
22.1520 
22.1577 

23.07 
25.00 

27.718 

35.22 



4.006 

2.665 

2.000 



1.6959 
1.6728 

1.6615 
1.6504 
1.6285 

1.335 
1.145 

1.001 

0.800 



**]Tor = 1 we used in part Southwell & Skan's data. 
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to equation 
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i 
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0.30 
0.70 
0.80 
0.86 
0.90 
0.S5 
0.99 
1. 01 
1.03 
1.50 
1.66 
2.00 



0o85? 
0.90 
0,93 
0.96 
0.S9 
1.05 
1.11 
1.20 



0.988 

0.9451 

0.9288 

0.9180 

0.9103 

0*9002 

0.8917 

0.8874 

0.8828 

0.7458 

0.6795 

0.4845 



1.810 

1.8956 
1.9138 
1.924 2 
1.9310 
1.9392 
i 1.S455 
j 1.9487 
l 1,9518 
'2.015 
2, ,032 
2.056 



26 



7.77 
5.365 
5.0596 
4.9014 



0.9101 
0.9037 
0.8970 



0.8899 
0.8827 
0.8674 



0.8509 
0.8237 



4.8039 
4.6930 
4. 6034 
4.5620 



4.5214 
3.7660 
3.5560 
3.1600 
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20 
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1,8462 
1.8484 
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1.8526 
1.8545 
1.8581 
1.8612 
1.8651 



4.5421 
4.4622 
4.3853 



0.524 
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0.4489 
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0.4344 

0.4272 

0.4213 

0.4185 

0.3406 

0.3097 
0.2201 



4.3114 
4.2397 
4,1034 



3.9751 
3.7948 
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0.4755 
0.4718 
0.4681 
0.4644 

0.4520 

0.4317 



20 



1.346 
1.4307 
1.4626 
1.4839 
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1.4992 
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1.5369 
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1.660 10.46 



14„08 
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10.823 

10.795 
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10.8741 
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13.390 
15.535 
28.89 



27 



4.265 
2.624 
2.402 
2.283 
2.209 
2.126 
2.054 
2.022 

1.361 
| 1.166 
I 0.735 
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12 
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Rigid Restraint 



20 



0.568 

0.5358 

0.5228 
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0.5024 
0,4964 
0 , 495 4 1 
0.4902 | 
0.4010 1 
0.3603 
0.249 
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10.07 



1.955 



1.810 



1.851 
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1.9504 

1.9843 

2.0035 

2.0199 

2.0284 

2.0375 

2.362 

2.568 

3.496 



0.50 62 



2.0026 



0.4916 | 2.0383 
0 .4811 I 2.0694 
2.1041 
2.1645 



C.4699 
0.4518 



21 



26.62 
19.77 
19.247 

18.915 

18.824 

18.783 

18.770 

18.771 

20.99 

23.45 

38.70 



17.895 



17.6141 

17.5716 



17.5991 
17.7792 
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3.096 
1.956 
1.801 

1.669 
1.608 
1.563 
1.540 



1.519 
1.084 
0.947 
0.635 
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TABLE I ( 



cont .) 
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2 


3 


4 


5 


6 1 


7 1 


8 


9 


10 I 


11 I 


12 


13 












Free support 


S-igid restraint 






n 


m 

m = — 
— i 
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i 


a 


Ca 


n i 

^a 


ts.. 

2. 


a 




Ca 


Ca' 


According 2 2 
to equation 


22 


26 


20 


20 


21 


27 


20 


20 


21 


27 




0.20 
0.30 
0.40 
0.50 
0.75 
0.92 
0.95 


0.995 

0.989 

0.980 

0.970 

0,93 

0.8926 

0.8850 


1.752 

1.766 

1.732 

1.779 

1.7790 

1.7785 


9.079 

7 442 

6.46 

5.78 

5.38 

4.1302 

4.0490 


0.5623 
0.5552 
0.5475 
0.54 

0.5048 
0.5018 


1.321 
1.332 
1.351 
1.365 

H Add C~> 

1 .495 8 
1.5100 


15.849 
13.205 
11.80 
10.744 

9 .2426 
9.2290 


5.317 
4.306 
3.679 
3.255 

2 . lyou 
2.1338 


0.5735 
0.5689 
0.5623 
0.5541 
0.5244 


1.840 

1.850 

1.8635 

1.8815 

1.9487 


30.747 
25.465 
22.415 
20.421 
17.69 


3.818 
3.101 
2.666 
2.361 
1.863 


5 


0.98 
1.00 
1.075 


0.8770 
0.8716 
0.8495 


1.7780 
1.7775 
1.7753 


3.9703 
3.9192 
3.7340 


0.4987 
0 . 49 65 


1.5256 
1 . 53 83 


9.2403 
9 .2743 


2.0801 
1 .959 6 


0 . 4841 
0.4691 


<3.0609 
2.1071 


16. 647 
16.586 


l,oe,5 
1.438 




X . XD 

1.25 
1.50 
1.75 
1.80 
1.85 


0.787 
0.666 
0.467 
0.408 
0.336 


1.7720 

1.766 

1.745 

1.708 

1,699 

1.690 


3.5621 

3.347 

2.S50 

2.371 

2.275 

2.180 


0.4635 
0.412 

0.2744 
0.2300 


1.726 
2.066 

3.5235 
4.335 


9.965 
12.155 
ci .y xy- 

28.248 
41.00 


1.628 
1 .23 8 
U. fol 
0.665 
0.533 


0.4523 
0.4277 
0.3521 
U.iooco 
0.2062 
0.1687 


2.1623 
2.2527 
2 . 6127 

T CI T a 
<J • DliO 

4.1011 
4.925 


16.655 

16.98 

19.43 

38.27 
52.9 


1.355 
1.248 
0.982 

U*DDC 

0.571 
0.469 




0.96 
0.S7 
0.98 


0.8759 
0.3730 
0.8701 


1.7203 
1.7194 
1.7184 


3 . 8032 
3.7755 

3.7477 


0.5321 
0.5312 
0.5305 


1.5113 
1.5168 
1.5222 


8.6865 
8.6861 
8.6836 


2 .1216 
2.1030 
2.0848 










10 


0.99 
1.00 
1.07 
1.08 
1.09 


6.8571 
0.8641 
0.8415 
0.8381 
0.8346 


1.7175 
1.7166 
1.7094 
1.7083 
1.7072 


3.7204 
3.6932 
3.5094 
3.4839 
3.4587 


0.5296 
0.5288 


1.5279 
1.5336 


8.6854 
8.6859 


2.0665 
2.0485 


0.4630 
0.4606 
0.4583 


2.1244 
2.1319 
2.1396 


2. 5 • 83 75 
15.8347 
15.8340 


1.4297 
1.4180 
1.4064 




1.10 
1.11 


0.8310 
0.8274 


1.7060 
1.7049 


3.4336 
3.4086 










0.4559 
0.4534 


2.1476 
2.1557 


15.8358 
15.8400 


1.3948 
1.3832 
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2 


3 


4 


5 


6 


7 


8 


9 


10 


11 
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n 
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straint 
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a 
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a 


C a 
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22 


22 


26 


20 


20 


21 


27 


20 


20 


21 


27 


20 


0.80 
0.95 
i.00 


0.914 
0.876 
0.86O 


1.708 
1. 682 
1.685 


4.1949 
3 . 7223 
3.5792 


0.5570 
0.5472 
0.5470 


1.438 
1.504 
1.533 


8,66 
8.43 
8.42 


2.44 

2.15 
2.05 


0.4751 

0.4562 
0.453.3 
0 . 4460 
0.4378 
0.4234 


2.084 

2.144 
2.160 
2.178 
2.207 
2.260 


15.549 

15.458 

15 . 448 
15.463 
15.499 
15.618 


1.507 

1.410 
1.386 
1.363 
1.328 
1.299 


1.05 
1.08 
1.10 
1.12 
1.15 
1.20 


0.843 

0.8325 

0.825 

0.817 

0.805 

0.785 


1.679 

1.674 

1.671 

1.667 

1.662 

1.'653 



3.4430 
3.3628 
3.3098 
3.2590 
3.1828 
3.0600 


0.5426 

0.5270 


1.566 
1. 690 


8.44 
8.71 


1.96 
1. 693 


40 


0,98 
1.00 


0.8643 
0.8578 


1 . 6721 
1.6691 


3.5790 
3.5217 


0.5572 
0.5560 


1.5199 
1,5320 


8 .,2678 
8,2655 


2.0879 
2.0507 


0.4514 
0.4539 
0.4487 
0.4433 
0,4349 


2.1256 
2.1501 

2,1675 
2.1854 
2*2154 

i 


15.2823 

15,2606 
15.2584 
15.2597 
15.2984 


1 . 4424 
1.4060 
1.3820 
1 .3584 
1.3224 


1.02 
1.05 

1.08 
1.10 
1.12 
1.15 


0.8511 
0.8407 
0.8297 
0.8220 
0 . 8141 
0.8016 


1 . 6660 
1.6611 
1.6560 
1.6525 
1 . 6488 
1.6432 


3 . 4653 
3.3824 
3.3011 
3.2478 
3.1951 
3.1170 


0.5546 


1.5447 


8,2682 


2.0138 
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TABLE I (cont.) 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


! 11 


12 


13 














Free 


Support 




Bigid H 


estraint 






Sa 


n 


m 

m = — 
_ 1 


* a 


_ % 

i 


a 


C 

a 


a' 


i 


a 


C a 


a' 


According 
to equatio 


22 

n 


22 


26 


20 


20 


21 


27 


20 


20 


21 


27 




0.15 
0.25 
0.50 
0.75 
0.95 
0.98 
1.00 


0.996 
0.992 
0.9675 

0.872 

0.8623 

0.8556 


1.730 
1.727 
1.714 

J. . Doi? 

1.661 

1.6563 

lo6527 


10.2918 
S.9720 
5.4758 

3.6114 
3.5224 
3.464 


0.577 
0.577 

0.577 

KJ . u / *x 

0.5684 
0.5667 
0.5656 


1.320 
1.319 
1.353 

1 AOl 

1.509 
1.51S 

1,5315 


17.96 

13.6 

10.00 

O. Co 

8.139 
8.128 
8.125 


6.15 

4.7625 
3.295 

Ci .O (K) 

2.147 

2.0896 

2.0515 


0.5743 
0.5707 
0.56 
0.52 

0.4714 


1.839 
1.843 

.1.887 
1.96 

2,092 


34.82 
27.00 
19.52 

1 R AA 

ID. '±'± 

15.14 


4.4075 

3.405 

2.36 

1 OC;7 
1 . ODO 

1.5025 


v. <■? 


1.02 
1.05 
1.08 
1 . 10 


0.8487 
0.8379 
0.82 66 
0.8185 


1.6493 
1.6438 
1.6381 
1.6340 


3.4059 
3.3222 
3.2395 
3.1846 


0.5644 
0.5625 


1.5445 
1.5655 


8.127 
8.142 


2.014 
1.958 


0.4591 
0.4514 
0.4458 


2.1312 

2.1566 
2.1750 


15.089 

15.067 
15.065 


1.4385 
1.4017 
1.377 


i 

1 


1.12 
1.15 
1.25 
1.40 
1.50 
1.60 
1.65 
1.66 
1.70 


0.8105 

0.7975 

0.749 

0.654 

0.568 

0.447 

0.359 

0.342 

0.333 i 


1.630 

1.624 

1.600 

1.560 

1.524 

1.483 

1.468 

1.453 I 

1.433 


3.1312 
3.0517 
2.7924 
2.4138 
2.1604 
1.8974 
1.7579 
1.7339 
1.6168 


0.540 

0.505 

0.462 

0.3866 

0.3218 

0.3062 

0.214 j 


1.768 

2.045 

2,388 

3.100 

3.945 

4.16 

6.32 j 


8.74 
10.10 
12.33 
18.23 
27.3 
29.95 
64. 65 


1.573 

1.299 

1.073 

0.801 

0.62 

0.586 

0*382 j 


0.4403 
0.4317 
0.3993 
0.3397 
0.2887 
0.2213 
0.1758 


2.1936 

2 • 225 

2.352 

2.652 

3.015 

3.76 

4.62 


15.067 

15.108 

15.43 

16.95 

19.68 

26.73 

37.35 


1.3535 

1.3165 

1.195 

1.000 

0.85 

0.66 

0.53 
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TABLE II 

Recapitulation of shear strength factors: c a and c\, (for 
computing t-^- r according to la and l"b) and c a ' and c-^' 

(for computing the confcrmpl wave length l^. T of the 
elastic surface according to 4a and 4b) . 
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2.93 
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19.9 
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1.36 
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2.49 


22.15 


22.15 


1.66 


1.66 


2 
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15.25 


2.28 


3.23 
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26.55 


1.54 


2.18 
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9.95 
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3.75 
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1.48 


2.57 
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Figs. 1,3a, 2b 
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Fig.l Elastic 
surface 
(contours and 
sections) of 
orthogonal 
anisotropic plate 
strip when wrinkling. Def lection w varies periodically in the 
longitudinal direction of the strip. The distance of the 
ensuing junction lines(the half-wave length) is 1. 
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Fig. 2a c a plotted 
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Figs.3a,3lD 
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Fig. 3b Quotient c^/c^id- = 0) plotted against l/£. 
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Figs. 4 a , 4b 
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Figs . 5a, 5b 
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Fig. 5a c a plotted against 1/& within range of 0 = l/S- = 1 
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Figs. 6a, 6b, 7 



m, 1.4*t 
n,2.13 

g, 2.49 

h, 1.66 

i, 2.28 
j,1.54 
k,1.48 
1,2.16 

Fig. 6a c a ' 

plotted 

against 



Curve 
ii 



'a 



\~g 



■?.M 



01 



El 1 



-6 



e, rigid restraint 
suoport 
3 



f , free 



e-i 



il 



10 20 
Fig. 6a 



in 

o 



"A 



03 



£>-g 




r 


i — 


w 














1.167 





30 0 



2 5 



10 
l/£ 



Fig. 6b 



g, 2.49 

h, 1.66 
.07 



1 « ^ 



j.1.36 
k,1.94 
1,1.20 



Fig. 6b C-J3 1 
plotted 
against l/£ 



20 



t 



HI 


ID 




FIT 








! 


! 


lo 
Si 




1 

> 




1 




— 


i 
1 


| 1 






| 

__j — 




-ri 

1 


A 


1 






1 

J t 








LO ( 
















! ll 

1 1 


j | 










t 


II C\2 

t\ 




1 


d=4 


r t ] 




6 = 0 
i-t- j 



X 



. 5 cm 



Fig. 7 Section of reinforced plate (example). 
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Figs .8a, 8b 
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Figs. 9a, 9b 
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Fig. 9b XiQ=(ri 3 -ra 2 ) fft,^ plotted against £ b and£,(0 = £ = l). 
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Figs- 10,11 
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Figs. 13a, '13b 
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Fig. 12a tyf/i and ^ e /i(>!/ e /l 2 ) plotted against £ a and 
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Figs. 13, 1-1 
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Figs . 15a, 15b 
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